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ABSTRACT 

We assess the practicality of computing the distance to stellar streams in ou r Galaxy, 
using the method of Galactic parallax suggested by lEvre fc BinnevI (|2009f ) . We find 
that the uncertainty in Galactic parallax is dependent upon the specific geometry of 
the problem in question. In the case of the tidal stream GD-1, the problem geometry 
indicates that available proper motion data, with individual accuracy ~ 4masyr~^, 
should allow estimation of its distance with about 50 percent uncertainty. Proper 
motions accurate to ~ Imasyr"^, which are expected from the forthcoming Pan- 
STARRS PS-1 survey, will allow estimation of its distance to about 10 percent uncer- 
tainty. Proper motions from the future LSST and Gaia projects will be more accurate 
still, and will allow the parallax for a stream 30 kpc distant to be measured with ~ 14 
percent uncertainty. 

We demonstrate the feasibility of the method and show that our uncertainty esti- 
mates are accurate by computing Galactic parallax using simulated data for the GD-1 
stream. We als o appl y the method to actual data for the GD-1 stream, published by 
iKoposov et al.l (|2009t ). With the exception of one datum, the distances estimated us- 
ing Galactic parallax match photometric estimates with less than 1 kpc discrepancy. 
The scatter in the distances recovered using Galact ic parallax is v e ry low , suggest- 
ing that the proper motion uncertainty reported by iKoposov et al.l (|2009t ) is in fact 
over-estimated. 

We conclude that the GD-1 stream is (8 ± 1) kpc distant, on a retrograde orbit 
inclined 37° to the plane, and that the visible portion of the stream is likely to be 
near pericentre. 

Key words: astrometry - methods: numerical - methods: data analysis - Galaxy: 
structure 



; 1 INTRODUCTION 

, Measuring distances in our Galaxy is critical to understand- 
ing its structure. However, line-of-sight distances can typi- 
cally be measured with only relatively poor precision. This 
lack of precision is manifest in the most basic of Galactic 
parameters; the solar radius Ro is hardly known to bet- 
ter than 5 percent uncertainty (|Gillessen et al.ll2009l l. and 
this res ult renders the circular vel ocity at Rq similarly un- 
certain (|McMillan fc Binnevll2009l ). An accurate knowledge 
of distances is essential to create convincing models of the 
Milky Way, which in turn influence our understanding of the 
physics of galaxy formation in general. 

Conventional trigonometric parallax has long been used 
to calculate accurate distances to nearby stars. The regular 
nature of the parallactic motion of a star, caused by the 
Earth's orbit around the Sun, allows this motion to be de- 
coupled from the intrinsic proper motion of the star in the 



heliocentric rest frame. Hence the distance to the star can be 
calculated. However, the maximum baseline generating such 
parallaxes is obviously limited to 2AU. For a given level 
of astrometric precision, this imposes a fundamental limit 
to the observable distance. Indeed, the accu racy of paral- 
laxes reported by the Hipparcos mission data (|van LeeuwerJ 
[2007) falls to 20-30 percent at best for distances ~ 300 pc 
and only then for the brightest st ars. Upcoming astr ometric 



projects such as Pan-STARRS feaiser et al.l I2OO2I') LSST 
(iTysonllioO^ ) and the Gaia mission (jPerrvman et al]l200ll ) 
will achieve similar uncertainty for Sun-like stars as distant 
as a few kpc, and at fainter magnitudes than was possible 
with Hipparcos. This extended range will encompass less 
than 1 percent of the total number of such stars in our 
Galaxy. 

It is clear that it will not soon be possible to calculate 
distances to many of the stars in our Galaxy with conven- 
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tional trigonometric parallaxes. Alternative means to com- 
pute distances to stars are therefore required. Photometry 
can be used to estimate the absolute magnitude of a star 
which, when combined with its observed magnitude, allows 
its distance to be computed. Unfortunately, all attempts to 
calculate such photometric distances are hindered by the 
same problems: obscuratio n by intervening m atter alters 
both observed magnitude (IVergelv et al. 1998) and colour 
l|Schlegel et al.lll998l : lDrimmel fc Soergel 20011 and it is dif- 
ficult to model appropriate corrections without a reference 
distance scale. The effects of chemical comp osition and age 
further complicate matters (|juric et al.ll200^ . It is therefore 
difficult to compute photometric distances with an accuracy 
much better than 20 percent, even for nea rby stars, and dis - 
tances to faint stars are less accurate still Jjuric et al]l2008h . 

Latterly, it was realised (|Evre fc Binnevll2009l . hereafter 
Paper I ) that the orbital motion of the Sun about the 
Galaxy could be used to compute trigonometric distances to 
stars. In the general case, it is not possible to do this because 
the parallactic motion of a star and its intrinsic proper mo- 
tion are inextricably mixed up. However, in the special case 
where the star can be associated with a stellar stream, its 
rest-frame trajectory can be predicted from the locations of 
the other associated stars. Using this trajectory, the proper 
motion in the Galactic rest frame can indeed be decoupled 
from the reflex motion of the Sun, and the component of its 
motion due to parallax can be computed. 

Such 'Galactic parallaxes' have the same geometrical 
basis as conventional trigonometric parallaxes, and as such 
as free from errors induced by obscuration and reddening. 
However, the range of Galactic parallax significantly exceeds 
that of conventional parallax. This is because the Sun orbits 
about the Galactic centre much faster than the Earth orbits 
the Sun, and because, unlike with conventional parallax, the 
Galactic parallax effect is cumulative with continued obser- 
vation. In realistic cases, for a typically oriented stream, we 
can expect the Galactic parallax to be observable at nearly 
40 times the distance of the equivalent trigonometric paral- 
lax, based on 3 years of observations. For increased range, 
one simply observes over a longer baseline. 

This large range means that Galactic parallax might 
prove a powerful tool to complement conventional paral- 
laxes, and validate other distance measuring tools. It is ex- 
citing to note that the capabilities of astrometric projects 
such as LSST, which will observe the conventional parallax 
of a G star at a distance of ~ 1 kpc with 20 percent un- 
certainty, will put much of the Galaxy in range of Galactic 
parallax calculations with similar accuracy. 



The main restriction on the use of Galactic paral- 
lax is the requirement for stars to be part of a stream. 
However, t he continuing discovery of significant numbers 



of streams todcnkirchcn et al. 2002|: iMaiewski et al.| 
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tical surveys implies that they are a staple feature of the 
Galactic environment, rather than a rarity. The deep sur- 
veys of Pan-STARRS and LSST are likely to find yet more, 
increasing the number of applications for Galactic parallax. 

This paper explores the viability of using Galac- 
tic parallax to estimate distances and demonstrates its 
practicality by applying it to data for the GD-1 stream 



(I Gr illmair fc DionatosI |2006| ) published by IKoposov et al.l 
(200l' hereafter K09 ); we choose to work w ith the latter 
over the earlier analysis of the same stream bv lWillett et al.l 
( 2009) on account of the significantly smaller proper motion 
uncertainties (Imasyr"^ vs 4masyr~^) cited in the later 
work. Throughout this paper, the Solar motion is assumed 
to be {U,V, W) = (10.0,252,7 1) ± (0 3, 11, 0.34) kms-\ 
consistent with Aumer & Binnov (l 2009l ). lReid fc Brunthalen 
(2004 ') and .Gillessen et al. (2009) - 

The paper is arranged as follows. Section [2] reviews the 
calculation of Galactic parallax, and explores the uncer- 
tainty affecting such calculations, and how this uncertainty 
affects practical application. Section [3] demonstrates the vi- 
ability of the method by applying it to pseudo-data, and 
Section [4] applies the method to actual data for the GD-1 
stream. Section [5] summarizes our conclusions. 



2 GALACTIC PARALLAX 

Suppose that a star is part of a stellar stream, and has a lo- 
cation relative to the Sun described by (x — xo) = rr, where 
r is the distance to the star, and xo is the position of the 
Sun. In the plane of the sky, let the tangent to the trajec- 
tory of the stream, near the star, be indicated by the vector 
p. Assume the velocity of the Sun, vo, in the Galactic rest 
frame (grf) is known. Paper I showed that if the measured 
proper motion of the star is /it, then 

tip = /it -I- — = /it + flVj,, (1) 
r 

where 11 = 1/r is the Galactic parallax, ii is the proper 
motion as would be seen from the grf, and Vs is the Sun's 
velocity projected into the plane of the sky. We note that u — 
vt/r, where vt is that component of the star's grf velocity 
perpendicular to the line-of-sight, and that, 

Vs = (vo ~ f • Vo f). (2) 

Eq.[T]is a vector expression and can be solved simultaneously 
for both u and H provided that p, t and Vj, are not parallel. 
The stream direction p will not typically be known outright, 
but must be estimated from the positions of stream stars 
on the sky. We can achieve this by fitting a low order curve 
through the position data, the tangent of which is then taken 
to be p. The curve must be chosen to reproduce the gross 
behaviour of the stream, but we must avoid fitting high- 
frequency noise, because p is a function of the derivative of 
this curve, which is sensitive to such noise. 

2.1 Uncertainty in Galactic parallax calculations 

We begin by rendering Eq. [T] into an orthogonal on-sky co- 
ordinate system, whose components are denoted by {x,y). 
In this coordinate system Eq. [T]can be solved for If, 

„ It (tx sin a — ty cos a) , . 

n = ^-^ ^ '-, (3) 

Vs^y COS a — Vs,x sm a 

where the (a;, y) suffixes denote the corresponding compo- 
nents of their respective vectors, and where we have defined 
the angle a = arctan(py/pi). 

The choice of coordinates {x,y) is arbitrary. We are 
therefore free to choose the coordinate system in which 
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a = 0, i.e. that system in which the a;-axis points along 
the stream trajectory, p. Eq. [3] becomes, 

n = -^, 

where we now identify the j/-component of the various vec- 
tors as that component perpendicular (±) to the stream 
trajectory, and the a::-component as that component parallel 
(II) to the trajectory. Eq.U shows explicitly that the Galac- 
tic parallax effect is due to the reflex motion of stream stars 
perpendicular to the direction of their travel. 

Uncertainties in the (a;, y) components of the measured 
quantities and Va, and uncertainty in a, can be propa- 
gated to n using Eq. [3] When we set a = 0, this equation 
becomes. 



+ ■ 



+ 



(5) 



where we anticipate the uncertainty in fit to be isotropic, 
and so we have set cr^jt^ — a^j^ty — a^. 

We assume to be known from observations; it may 
contain any combination of random and systematic error. 
av^j^ is calculated directly from the error ellipsoid on vo, 
which is assumed known. Any error on vo affects all data in 
exactly the same way. However, the projection of error on 
Vo to Vs varies with position on the sky. Hence, the effect 
of CTv^ is to produce a systematic error in reported distance 
that varies along the stream in a problem-specific way. 

Uncertainty in a arises from two sources. Firstly, be- 
cause the on-sky trajectory p is chosen by fitting a smooth 
curve through observational fields, p need not be exactly 
parallel to the underlying stream. Further, since p depends 
on the derivative of the fitted curve, it is likely to be much 
less well constrained for the data points at the ends of the 
stream than for those near the middle. 

We can quantify this effect. At the endpoints, the fitted 
curve is likely to depart from the stream by at most A?/), 
the angular width of the stream on the sky. For a low-order 
curve, this departure is likely to have been gradual over ap- 
proximately half the angular stream length, AO, giving a 
contribution to from fitting of. 



2 



(6) 



The second contribution to aa arises as follows. Since 
the stream has finite width, at any point, the stars within it 
have a spread of velocities, corresponding to the spread in 
action of the orbits that make up the stream. If the stars in 
a stream show a spread in velocity {av^ , (Jvy ) about a mean 
velocity vt = riip, this effect contributes. 



2 1/2 2 2.2\ 

^a.v = ~ Wv COS OL -f (Tj,^ Sm Q j 



(7) 



to the uncertainty in a for a single star. Again we can choose 
Q = 0, such that a^y ~ o-v±, the velocity dispersion perpen- 
dicular to the stream direction. Eq. [7] becomes. 



(8) 



(iisin/3)2' 

where we have introduced v, the grf speed of the stream, 
and /3, the angle of the stream to the line-of-sight. av± has 
its origin in the random motions of stars that existed within 
the progenitor object. In fact, if we assume the stream has 



not spread signific antly in width, then the w idth and the ve- 
locity dispersion (|Binnev fc Tremainell2008l . §8.3.3) are ap- 
proximately related by. 



w 



rAtp 

Rn 



(9) 



where w is the physical width of the stream, and Rp is the 
radius of the stream's perigalacticon. This gives, 

rAtp 



Rp sin /3 



(10) 



If secular spread has made the stream become wider over 
time, then this relation will over-estimate aa^v, since a„j_/v 
is roughly constant. Ail) therefore represents an upper bound 
on the true value of aa.v through this relation. This argu- 
ment also assumes that the stream was created from its pro- 
genitor in a single tidal event. Real streams do not form in 
this way. However, repeated tidal disruptions can be viewed 
as a superposition of ever younger streams, created from a 
progenitor of ever smaller a^x. Eq. [10] holds for each of these 
individually. Thus, Atp remains a good upper bound for aa,v 
through this relation. 

In reality, we do not measure the proper motion of in- 
dividual stream stars, but rather the mean motion of a field 
of A*' stars. The contribution to is from the error on this 
mean. Putting this together with Eq. [6] gives our final ex- 
pression for CTq, 



N 



+ ^c,f = 



-I- 



AAip 



NRl sin'^P AO 



(11) 



We note that the first term represents a random error, and 
the second term represents a systematic error that will vary 
with position down the stream. In general, sin (3 and _Rp are 
a priori unknown. We can infer sin/3 from radial velocity 
information, either directly where the measurements exist, 
or indirectly from Galactic parallax distances. Guessing Rp 
requires assumptions to be made about the dynamics, but 
in general we expect the ratio r/Rp ~ 1 or less. 

Explicit evaluation of sin 13 and Rp are not necessary 
to evaluate the uncertainty if ctq is dominated by the error 
from fitting, ctq,/. We can see this will be the case when the 
number of observed stars per field. 

We expect this to be true in almost all practical cases. 

2.2 Uncertainty in tangential velocity calculations 
Eq. [T]can also be used to solve for ii, 

fl{ty +tx) +'n-{Vs,y + Vs,x) 

u = ; — ■■ , 

cos a + sm a 

which becomes, 



(13) 



il = /i(t|| + t±) + U{Vs\\ + Ua±) = + UVs 



(14) 



when we set q = 0. Eq. [13] combined with Eq. |3] can be 
used to explicitly propagate uncertainties in the measured 
quantities to ii. When a = 0, the uncertainty in ii is, 

aj _ ^ ^i(«.\<|| +"f||<^) 



© 2009 RAS, MNRAS 000, [TH9] 



4 A. Eyre 



2t]_Vs\\Vs±COY{Vs\\,Vs±) , VsW^l , , 

— 9 7 \o 9 ■ 

cr„_,|| and cov(?;s|| , Usx) are calculated directly from the error 
ellipsoid on vo, which we have assumed known. 

2.3 Practicality of Galactic parallax as a distance 
measuring tool 

Using Eq.Uto eliminate fit± from Eq.[Sl and taking the dot 
product of p with Eq. [T]to simplify the last term, we obtain, 

'^n 1 fc \2 , 2 , / .n2 2I 

= + (16) 

where we have noted that rii — Vt = v sin /3, and we have 
related the observed stream length, A6, to the deprojected 
length, A9 = AS sin /3. We note that the last term is inde- 
pendent of r, since rAip = w and Aip/AQ are both constant, 
and that for a stream of given physical dimension, the uncer- 
tainty in n has no dependence upon the angle of the stream 
j3 to the line of sight. 

What level of uncertainty does Eq. [16] predict, when 
realistic measurement errors are introduced? The answer to 
this is dependent upon the both the physical properties of 
the stream {Rp, Atp, AO, v) and the geometry of the problem 
in question {r,Vs±). 

We progress by assuming 'typical' values for some of 
these quantities. The average magnitude of taken over 
the whole sky is vqtt/4. The average perpendicular compo- 
nent, for a randomly oriented stream, is 2/tt of this value. 
We therefore assume a typical value for Vs± of uo/2 
120kms"\ We also assume a typical grf velocity equal to 
the circular velocity, v = Vc ~ 2 20kms~^. 

iMcMillan fc BinnevI l|2009l 'l recently summarised the 
current state of knowledge of vq. The uncertainty quoted 
is typically --^5 percent on each of {U,V,W). Correspond- 
ingly, we estimate a typical value for the uncertainty (J„^j^ 
of 5 percent of or 6kms~^. 

The GD-1 stream that we consider below is exception- 
ally thin and l ong, with At/; ~ 0.1° and A9 ~ 60° . The 
Orphan stream (|Grmmaij|2006l : iBelokurov et al.ll2007^ is of 
similar length, but about 10 times thicker. Both of these 
streams are near apsis, so A9 = AO. We therefore take 
Alp 1°, AO ~ 60° as typical of the streams to which one 
would apply this method. If hundreds of stars are observed 
for each proper motion datum, then Eq. [12] is true for all 
realistic combinations of (r, 7?p), so we can ignore the con- 
tribution of aa,v to cTa. The contribution from a^j gives 
(j<,~1.9°. 

The in dividual USNO/SDSS proper motions 
l|Munn et al.l 2004) used by K09 have a random un- 
certainty CTp ~ 4masyr~^. After averaging over hundreds 
of stars and accounting for a contribution from non-stream 
stars, K09 report a random uncertainty of ct^ ~ Imasyr"^ 
on their GD-1 data. For a stream lOkpc distant, with 
these proper motions and the typical values mentioned, 
Eq. [TS] reports an uncertainty of an/II ~ 40 percent. By 
far the greatest contribution comes from the first term in 
Eq. 1161 hence, the error on proper motion measurement is 
dominating our uncertainty. 



To obta i n an uncertainty of crn/II < 20 percent with 
iMunn et all (|2004l ) proper motion measurements, we would 
need to restrict ourselves to streams less than 5kpc dis- 
tant. 20 percent error is also possible at lOkpc given opti- 
mum problem geometry. This is clearly competitive with the 
^ 20 pc at which one could observe a standard trigonomet- 
ric parallax, with similar accuracy, using astrometry o f this 
quality. However, previous work (^W illett et al.l l200Sj -. K09 
) shows that SDSS photometry combined with population 
models produce distance estimates accurate to ~ 10 per- 
cent for stars in streams at 8kpc. The accuracy of Galactic 
parallax is therefore not likely to be as good as that of pho- 
tometric distances for distant streams, using data this poor, 
unless the problem geometry is favourable. 

Proper-motion data from the Pan-STARRS telescope 
is expected to be accurate to ~ l masyr~^ for Sun- like 
stars at lOkpc (|Magnier et al.|[2008l) . K09 reduce raw data 
with accuracy ~ 4masyr~^ to processed data accurate to 
~ Imasyr"^, even though the expected proper motion of 
the stars is of the same size as the errors. It is not unreason- 
able to expect a similar analysis applied to Pan-STARRS 
raw data, where the relative error would be much less than 
unity, to yield processed data accurate to ~ 0.2masyr~^. In 
truth, the ability of Pan-STARRS to detect very faint stars 
will increase the number of stars identifiable with a stream, 
and thus reduce the uncertainty in the mean proper motion 
further than this, but we use 0.2masyr~^ as a conservative 
estimate. 

The same lOkpc distant stream would have a parallax 
error of crn/H ~ 11 percent with data this accurate. An er- 
ror of less than 20 percent is possible for a typical stream 
less than ~ 23 kpc distant, and for a stream w ith favourable 
geometry less than 50 kpc distant. I Juric et"al 

] I2OO8 ) report 

that SDSS photometric distances for dwarf stars have ~40 
percent error at 20 kpc. Thus, the accuracy of Galactic par- 
allax derived from Pan-STARRS data should be at least 
comparable to distance estimates from photometric meth- 
ods, even in the typical case. 

Future projects such as LSST and Gaia will each ob- 
tain proper motions acc urate to ^ 0.2masyr~^ for Sun- 
like s tars 10 kpc distant (|lvezic et al.ll2007l : jPerrvman et al.l 
I2OOII ). These data would allow a distance estimate for our 
typical stream accurate to 8 percent, and a stream with 
favourable geometry accurate to 4 percent. Error in the 
proper motion no longer dominates the uncertainty in these 
calculations. We might expect such accurate astrometric sur- 
veys to reduce the uncertainty in the Solar motion; in this 
case, the error in parallax would be lower still. 

Gaia will not observe Sun-like stars beyond 10 kpc, but 
LSST will, wi th accuracy of .4masyr~^ for dwarf stars 
30 kpc distant (jlvezic et al .I2OO7I ). The accuracy of the paral- 
lax to our typical stream at this distance would be about 14 
percent with these data, and 6 percent is achievable with op- 
timum geometry. A typical stream could be measured to 20 
percent accuracy out to 40 kpc, and a stream with favourable 
geometry out to 54 kpc; this range approaches the limit of 
LSST's capability for detection of dwarf stars. Such data 
will put the Orphan stream, which is about 20 — 30 kpc 
distant (iGriUmaiij l2006l : iBelokurov erall l2007l : ISales et al.l 
l2008h . in range of accurate trigonometric distance estima- 
tion. For comparison, photometric distances from SDSS data 
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Figure 1. Full line: the orbit for the GD-1 stream, taken from 
K09. Crosses: pseudo-data derived from that orbit, but randomly 
scattered in (l>2 according to a Gaussian distribution with a dis- 
persion of tT^2 = 0.1°. Dotted line: a cubic polynomial fitted to 
the pseudo-data, used to estimate stream direction. 
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Figure 2. Dotted line: the orbit of the GD-1 stream, taken from 
K09. Crosses: the true distance of each pseudo-datum. Circles: 
Galactic parallax distances computed from the pseudo-data. The 
error bars represent the distance error expected from the poly- 
nomial fitting procedure. No extraneous error was added. The 
estimated errors are show to be a good estimate of likely error 
from the fitting procedure, and the agreement of the distances 
overall is excellent. 



are hardly more accur ate than 50 percent for this stream 
l|Belokurov et al.ll2007^ . 



3 TESTS 

To test the method, pseudo-data was prepared from an orbit 
fitted to data for the GD-1 stream by K09. The orbit is de- 
scribed by the initial conditions x = (—3.41, 13.00, 9.58) kpc, 
V = (—200.4, —162.6, 13.9) kms~^, where the x-axis points 
towards the Galactic centre, and the y-axXs points in the di- 
rection of Galactic rotation. The orbit was integrated in the 
logarithmic potential, 



$(a;,j/. 




(17) 
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where Vc — 220 km s ^ and q = 0.9. The resulting trajec- 
tory was projected onto the sky, assuming a Solar radius 



Figure 3. Crosses: Galactic parallax distances computed from 
the pseudo-data, with no extraneous errors. Error bars: the ran- 
dom scatter expected in Galactic parallax distances, with mea- 
surement errors as mentioned in the text. Plus signs: Galactic 
parallax distances computed from 60 Monte Carlo realisations of 
each pseudo-datum convolved with the measurement errors. The 
analytic uncertainty estimate and the Monte Carlo realisations 
are in good agreement. 



Ro ~ 8.5 kpc. Several points were sampled, and each was 
taken to be a separate datum in the pseudo-data set. The 
proper motion for each datum was computed by projecting 
the difference between its grf motion and the Solar motion 
on to the sky. 

The pseudo-data were transformed into the rotated co- 
ordinate system used by K09 to facilitate comparison with 
their data; the transformation rule is given in the appendix 
to K09. The stream is very flat in this coordinate system, 
so the dependence of (j}2 on (j}i is relatively weak. This helps 
to increase the quality of the fitted curve and minimises the 
corresponding error in CTq,/. 

To simulate the observed scatter in the real positional 
data, the pseudo-data were each scattered in the (f>2 coor- 
dinate according to a randomly-sampled Gaussian distribu- 
tion with a dispersion ct^j — 0.1°. The resulting positional 
pseudo-data are plotted in Fig. [TJ along with the orbit from 
which they were derived (full curve). A cubic polynomial 
representing 4>2{4'i) was least-squares fitted to the pseudo- 
data, the tangent of which was used to estimate p. In the 
case of the pseudo-data, uniform weights were applied to 
each datum for the fitting processes. The resulting curve is 
also show in Fig. [T] (dotted curve). 

When the correct orbit is used to calculate p, and pre- 
cise values for the measured proper motion fit and Solar re- 
fiex motion Vs are used, the distance is recovered perfectly 
from Eq. U) Fig. [2] compares the recovered distance when 
p is estimated using the polynomial fit to the pseudo-data, 
but still using accurate values for fit and v^. Our pseudo- 
data stream is Aip ^ 0.1° wide and A9 ~ 60° long. Eq. [S] 
therefore estimates a^j — 0.38°. The recovered distances in 
Fig. [2] are in error by only ~ 2 percent across most of the 
range, which is the approximate uncertainty predicted by 
Eq. Ofor this value of iJaj- Thus, the estimation of p from 
the observed stream is good, and contributes little error to 
the distance calculations. 

The K09 observational data for the GD-1 stream, dis- 
cussed below, have a similar uncertainty ~ 0.38° due 
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Figure 4. Crosses: on-sky position data for tlie GD-1 stream, 
as published in K09. The error bars represent the quoted uncer- 
tainties. Pull line: linear least-squares fit of a cubic polynomial, 
<t>2{<j>i), to these data; the inverse-square of the uncertainties was 
used to weight the fit. 
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Figure 5. Circles: Galactic parallax distances for the GD-1 data 
presented in K09. Dotted error bars: the uncertainty estimated by 
Eq.O given the K09 measurement uncertainties. Crosses: the pho- 
tometric distances reported in K09, along with their error bars. 
Full line: the orbit for GD-1 taken from K09. With the exception 
of the datum at (f>\ ~ —55 deg, the Galactic parallax distances 
are in excellent agreement with the photometric distances from 
K09. The dotted error bars appear to seriously over-estimate the 
true error in the distance estimates. 



entirely to the fitting process, and proper motion uncertain- 
ties cr^ ~ Imasyr"^. Fig. |3] shows the recovered distances 
from Fig. [2] with error bars for the expected uncertainty in 
recovered distance, given these measurement uncertainties 
and the uncertainty in vo quoted in Section 1. Also plotted 
for each datum are the distances recovered from 60 Monte 
Carlo realisations of the pseudo-data input values, convolved 
with the errors given above. 

Eq.[5]is found to be a good estimator for the uncertainty, 
with approximately 80 percent of the Monte Carlo realisa- 
tions falling within the error bars. The error in parallax for 
the K09 data is thus predicted to be about 50 percent, of 
which the greatest contribution comes from the uncertainty 
in proper motion. 
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Figure 6. Full linos: Galactic-rest frame proper motion (ii) cal- 
culated from the K09 data using Eq. 1141 The (upper, lower) 
panels show the {(j>i,(j>2) components respectively. Plotted in the 
background are the observational data from Fig. 9 in K09; the 
greyscale shows the number of stream stars, per bin, with the 
given motion. The data are broadly consistent, except for the 
datum at cf>i ~ —55 deg in the upper panel. 

4 DISTANCE TO THE GD-1 STREAM 

Fig. m shows the on-sky position data for the GD-1 stream, 
as published in K09. Also shown in Fig. |4]is a linear least- 
squares fit of a cubic polynomial to these data, used to es- 
timate p. The weights for the fit were the inverse-square 
uncertainties for each position field, as given by K09. 

K09 provide measured proper motion data for five fields 
of stars, spanning the range 01 ~ (—55,-15)°, along with 
uncertainties for these measurements. Uncertainty in the 
stream direction is Oa ~ 0.38°, which is entirely contributed 
by the curve fit to the stream; since hundreds of stars con- 
tributed to the calculation of the proper motions, the con- 
tribution from the first term in Eq. [TT] is negligible. The 
uncertainty in Vs is computed for each individual field from 
the uncertainty in vo given in Section 1. 

Fig. [5] shows the Galactic parallax distances for each of 
these data, along with the K09 photometric distances. The 
dotted error bars represent the expected error in distance for 
the uncertainties given. The small solid error bars are the un- 
certainties reported by K09 for their photometric distances. 
The K09 orbit used to compute the earlier pseudo-data is 
plotted for comparison. 

With the exception of the datum at 01 ~ —55°, the 
parallax distances and the K09 distances are in remarkable 
agreement. However, the dotted error bars vastly overesti- 
mate the true error in the results. If we ignore the datum at 
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<j)i ~ —55°, the scatter in the distance, ar ^ 1 kpc, is similar 
to that of the photometric distances, and consistent with a 
true random error of cr^ ~ 0.3masyr~^, and neghgible sys- 
tematic offset. We cannot explain this discrepancy, except 
by suggesting that the K09 proper motion measurements 
are more accurate than the published uncertainties suggest. 
This is corroborated by the top-right panel of Fig. 13 from 
K09 in which the data, with the exception of the datum 
at (j>\ ~ —55°, show remarkably little scatter within their 
error bars. 

Fig. |S]shows the Galactic rest-frame proper motions, ii, 
calculated from Eg. ll4l along with their error bars, from Eq. 
1151 In the background are plotted the data from Fig. 9 of 
K09, which show the density of stars with a given grf proper 
motion in the sample of stars chosen to be candidate mem- 
bers of the stream, and after subtraction of a background 
field. The K09 grf proper motions have been calculated by 
correcting measured proper motion for the solar reflex mo- 
tion, using an assumed distance of 8 kpc (Koposov, private 
communication) ; this assumption will cause a systematic er- 
ror in the K09 proper motions, of order the distance error, 
which changes with position down the stream. The appar- 
ently large width of the stream in thi s plot is due to uncer- 
tainty in the underlying iMunn et al.l l|2004 ) proper motion 
data. 

The stream is clearly visible in this plot as the region of 
high density spanning <j)i ~ (0, —60)° with ~ Omasyr^^ 
and falling slowly between (—6, —10) masyr"^. Despite 
the expected systematic error, the estimates of ii from the 
parallax calculation are consistent with these data, with the 
exception of the same datum at cj>\ ~ —55° that also reports 
an anomalous distance. 

We explain this suspect datum as follows. From inspec- 
tion of the top-right panel of Fig. 13 from K09, it is apparent 
that the measurement for this datum is not in keeping 
with the trend. Conversely, the corresponding measure- 
ment is not obviously in error. If the magnitude of for 
this datum has been over-estimated by the K09 analysis, 
then Eq.|4]will over-estimate the parallax, and hence under- 
report the distance. Fig. [S] indicates that the distance for 
this datum is indeed under-reported. 

The effect of such an error in on the grf proper- 
motion, u, can be understood by considering Eq. 1141 If 11 
is over-estimated, ii will be either over-estimated or under- 
estimated, depending on the relative sign of the two terms. 
In the case of GD-1, /ity and Vs\\ have opposite signs, so an 
over-estimated II will result in an under-estimated ii. This 
too corresponds with the behaviour of the suspect datum in 
Fig. El 

It is unknown why this particular datum should be sig- 
nificantly in error while the other data are not. There are no 
obvious structures in the lower panel of Fig. |S] which might 
cause the fitting algorithm in K09 to mistakenly return an 
incorrect value for fi^^ . Nonetheless, if the scatter in the 
other data are accepted as indicative of their true statisti- 
cal error, it is clear that the datum at cj>\ ~ —55° cannot 
represent the proper motions of GD-1 stars at that location. 
We therefore predict that an appropriate re-analysis of the 
proper-motion data, taking care to ensure that a signal from 
GD-1 stream stars is properly detected, will return a revised 
proper-motion of ~ — 3masyr~^. 

In summary, it seems that Galactic parallax measure- 



ments confirm the K09 photometric analysis, and predict 
that the stream is approximately (8 ± 1) kpc distant, where 
the uncertainty denotes the scatter in the results. Since 
Galactic parallax and photometric estimates are fundamen- 
tally independent, it seems unlikely that systematic errors in 
either would conspire to produce the same shift in distance; 
this implies that no systematic error is present. 

We also calculate a grf proper motion for the stream of 
/^^j^ = ( — 7 ± 2) masyr"^, corresponding to a grf tangential 
velocity of (265 ± 80)kms~^ in a direction (/i; cos 6, /it) ~ 
(0.8, —0.6). This implies that the stream is on a retro- 
grade orbit, inclined to the Galactic pl ane by ~ 37°, which 
is in accordance wit h previous results l|Willett et al.ll2009l : 
iKoposov et"al]|2009l ). 

The galactocentric radius of ^ 14.5 kpc does not seem 
to be changing rapidly along stream's length, which sub- 
tends ~ 12° when viewed from the Galactic centre. This 
implies that the observed stream is at an apsis. The grf 
velocity of the stream is faster than the circular velocity, 
Vc ~ 220kms~^. This implies that the stream is at peri- 
centre, although the large uncertainty prevents a firm con- 
clusion from being drawn. We note that the radial velocity 
data in K09 would also imply that the stream is observed 
at pericentre. 



5 CONCLUSIONS 

We have demonstrated the practical application of a tech- 
nique for computing Galactic parallax, as described by Pa- 
per I. This technique utilises the predictable trajectories of 
stars in a stream to identify the contribution of the refiex 
motion of the Sun to the observed proper motion. The par- 
allax and the Galactic rest- frame proper motion follow from 
this. 



The only assumption made is knowledge of the Galac- 
tic rest-frame velocity of the Sun. It is also a require- 
ment that the observed stars are p art of a stream. Recent 



evidence ("Odonkirchon et al.l 


20021: 


iMaiewski et al' '2003; 


Yanny et al. 


2003; Bclokurov et al. 


[2006: Grillmair 2006; 


Grillmair & 


Dionatosi 120061: 


Grillmair & Johnson 12006: 


Grillmair |2009|; Newberg. Yanny, & WillettI 120091) indicates 



that tidal streams are a common constituent of the Galactic 
halo, and so this technique should have widespread applica- 
tion. 



We have derived an expression for the uncertainty 
in Galactic parallax calculations. We include contributions 
from measurement errors in proper motion and Solar mo- 
tion, error in the estimation of stellar trajectories from the 
stream direction, and algorithmic error in the estimation of 
stream direction itself. 

The uncertainty for calculations involving a particular 
stream is depend upon the size, location and orientation of 
the stream, as well as upon measurement errors. We estimate 
that using individual proper motions accurate to 4 mas yr~^ , 
available now in published surveys (|Munn et al.ll2004^ ■ the 
parallax of a 10 kpc distant stream with typical geometry 
could be measured with an uncertainty of 40 percent. The 
parallax of a stream with optimum geometry could be mea- 
sured with approximately half this uncertainty. 

Proper motion da ta from the forthcoming Pan- 
STARRS PS-1 survey (|Kaiser et al.l |2002| : iMagnier et"aLl 
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|2008| ) will yield the distance to a typical lOkpc distant 
stream with 11 percent accuracy, or the distance to a stream 
at 23 kpc with 20 percent accuracy; with favourable geome- 
try this accuracy could be achieved for a stream as distant 
as 50 kpc. With data of this quality, the uncertainty in dis- 
tances from Galactic parallaxes will be considerably lower 
than those of photometry for distant streams. 

The LSSrJ^son 2002': Jlvezic et al.ll2007l ') and the Gaia 
mission ijPerrvrnan et al..,20[)li ) will produce proper- motion 
data that are more accurate still. Such data would allow the 
distance to stars in a 10 kpc typical stream to be computed 
to an accuracy of order of 8 percent, where the limitation 
is now imposed by uncertainty in the solar motion and in 
the stream trajectory. It is likely that LSST and Gaia data 
will allow the uncertainty in the solar motion to be signif- 
icantly reduced, so in reality much better precision can be 
expected at this distance. For streams 30 kpc distant, LSST 
proper motions will allow distance estimates as accurate as 
14 percent to be made in the typical case, and 6 percent 
with optimum geometry. Thus, the high-quality astrometric 
data that is expected to be available in the next decade will 
allow parallax estimates for very distant streams to be made 
with unparalleled accuracy. 

To test the method presented, we have created pseu- 
doda ta simulating the GD-1 stream IjGrillmair fc Dionatod 
I2OO6I ). When the method is provided with error- free pseudo- 
data, the correct parallax is computed perfectly. When er- 
rors are introduced into the pseudo-data, the reported par- 
allax degrades in line with the uncertainty estimates. 

We applied the method to the astro metric data for 
the GD-1 stream in IKoposov etal\ (|2009l '). With the ex- 
ception of a single datum, the Galactic parallax is re- 
ma rkably consisten t with the photometric distances quoted 
by IKoposov et al.l (|2009l ). Indeed, t he uncertainty in the 
measured proper motions quoted by IKoposov et al.l (|2009| ) 
should produce significant error in the Galactic parallax. 
However, the scatter in the results is consistent a random 
error of only ~ 0.3masyr~^, and if the photometric dis- 
tances of K09 are believed, no systematic offset. This is at 
odds with the typical unc ertainty in the p roper motion of 
~ Imasyr"^ reported by IKoposov et all (|2009| ). We can- 
not explain thi s discrepancy, other than to suggest that the 
IKoposov et al.l (|2009h method for estimating error in the 
proper motions is producing significant over-estimates. 

The Galactic rest-frame proper motions predicted for 
the stream are also c onsistent with observational data from 
iKoposov et"al](|2009D . with the exception of the same datum 
that also reports an inconsistent distance. We conclude that 
the proper-motion associated with this datum is erroneous, 
and we predict that reanalysis of the stream stars near this 
datum will reveal a reduced proper-motion measurement of 

~ 3masyr~^. 

Photometry and Galactic parallax produce fundamen- 
tally independent es timates of distance. T he quality of the 
corroboration of the lKoposov et ahl (|2009l l photometric dis- 
tance estimates for GD-1 by the Galactic parallax estimates 
presented here therefore lends weight to the conclusion that 
the predicted distance, in both cases, is correct. On this ba- 
sis, we conclude that the GD-1 stream is about (8 ± 1) kpc 
distant from the Sun, on a retrograde orbit that is in- 
clined 37° to the Galactic plane with a rest-frame velocity of 



(265 ±75) kms~^. We also conclude that the visible portion 
of the stream is probably at pericentre. 

The prospect of being able to map trigonometric dis- 
tances in the Galaxy to high accuracy at tens of kiloparsecs 
range is indeed exciting. The distances generated using this 
method, although limited to stars in streams, could be used 
to calibrate other distance measuring tools, such as photom- 
etry, that would be more widely applicable. The technique 
is immediately applicable to any stream for which proper- 
motion data are currently available, although we anticipate 
limited accuracy until better proper-motion data are avail- 
able. 

Given enough parallax data points along a given stream, 
an orbit can be constructed by connecting those points. This 
orbit is predicted independently of any assumption about 
the Galactic potential, which it must strongly constrain. 
Constraints on the Galactic potential impose constraints on 
theories of galaxy formation and cosmology. It would seem 
that the combination of dynamics and Galaxy-scale preci- 
sion astrometry, such as provided by this method, could well 
have profound implications for astrophysics in the future. 
At present, however, it is not obvious how to combine all 
sources of astrometric and dynamical information, to pro- 
duce the tightest constraints on the potential. We therefore 
encourage the exploration of methods for combining this in- 
formation, in anticipation of the arrival of higher quality 
astrometric data in the next few years. 
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